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Abstract 

We consider the geometric random (GR) graph on the d— dimensional torus 
with the La distance measure (1 < a < oo). Our main result is an exact char- 
acterization of the probability that a particular labeled cycle exists in this random 
graph. For a = 2 and a = oo, we use this characterization to derive a series which 
evaluates to the cycle probability. We thus obtain an exact formula for the expected 
number of Hamilton cycles in the random graph (when a = oo and a — 2). We 
also consider the adjacency matrix of the random graph and derive a recurrence 
relation for the expected values of the elementary symmetric functions evaluated 
on the eigenvalues (and thus the determinant) of the adjacency matrix, and a re- 
currence relation for the expected value of the permanent of the adjacency matrix. 
The cycle probability features prominently in these recurrence relations. 

We calculate these quantities for geometric random graphs (in the a = 2 and 
a — oo case) with up to 20 vertices, and compare them with the corresponding 
quantities for the Erdos-Renyi (ER) random graph with the same edge probabili- 
ties. The calculations indicate that the threshold for rapid growth in the number of 
Hamilton cycles (as well as that for rapid growth in the permanent of the adjacency 
matrix) in the GR graph is lower than in the ER graph. However, as the number 
of vertices n increases, the difference between the GR and ER thresholds reduces, 
and in both cases, the threshold ~ log(n)/n. Also, we observe that the expected 
determinant can take very large values. This throws some light on the question of 
the maximal determinant of symmetric 0/1 matrices. 



1 Overview 

Consider the d— dimensional unit torus Td = [0, l] d . For < r < 1/2, 1 < a < oo, 
the geometric random (GR) graph Qn (r) is defined as follows. The vertex set cor- 
responds to n points X n = {xi, X2, ■ ■ - x n } distributed uniformly and independently 
in Td. The set of edges E{Q^' d \r)) is defined as 

E(Q { n' d) (r)) - {{x i ,x j }:\\x i -x j \\ q < r} 
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where || . || g is the L q norm. Then, Q ( y\r) is a random graph. In this random graph 
model, the presence of an edge is not necessarily independent of the presence of other 
edges. 

Another random graph model which has been very well studied is the Erdos-Renyi 
(ER) random graph, which is defined as follows. Given a number p, < p < 1, let 
H(n,p) denote the graph which has the vertex set {1, 2, . . . n} and an edge set con- 
sisting of edges selected with probability p (a particular edge {i,j} is present with 
probability p and the presence of each edge is independent of the presence of other 
edges). The ER random graph has been extensively studied. Specifically, the asymp- 
totic behaviour (or evolution) of this random graph has received considerable attention 
(DID. The most celebrated result of this type |1| can be summarized as follows: if 
p = p(n) = (log ri + c n )/n, then the random graph G„ is almost surely connected 
(as n — > oo) if c„ — > oo, and is almost surely disconnected if c„ — > — oo). Similar 
thresholds exist for all monotone graph propertied Q. 

The geometric random graph appears to exhibit similar asymptotic properties. In 
iflOl . a sharp threshold for connectivity has been exhibited for the geometric random 
graph on the unit square (d = 2 and a = 2): if r ~ r(n) and if 7rr(n) 2 = (logn + 
c n )/n then the random geometric graph is almost surely connected if c n — > oo, and 
is almost surely disconnected if c n —> — oo. The existence of sharp thresholds for 
monotone properties in geometric random graphs has been demonstrated in (51 . The 
monograph |5| summarizes threshold characterizations of several connectivity related 
properties of the geometric random graph. Upper and lower bounds on the diameter of 
a geometric random graph in the unit ball have been derived in (4). The mixing times 
of random walks in geometric random graphs have been characterized in O. The 
limiting distribution of the eigenvalues of the adjacency matrix of a random graph has 
been studied in Q, 0. An asymptotic bound for the second largest eigenvalue of the 
adjacency matrix of a geometric random graph has been derived in ifTTl . Thus, there is 
a large body of work on the asymptotic properties of a geometric random graph. 

In the finite case, one is interested in the exact formula for the appearance of a 
certain property in a geometric random graph. An example of such a characterization 
is an exact formula for the probability of connectivity of a geometric random graph 
on a 1-dimensional unit cube HI 211 . and an exact formula for the probability of exis- 
tence of a particular labeled subgraph in the geometric random graph constructed in 
the d— dimensional unit cube using the measure fl3l . We will consider the finite 
case, and prove an exact characterization of the probability that a labeled cycle appears 
in the random graph Q„ (r) (valid for 1 < a < oo, and for all d > 1). Using this 
characterization, we show that it is possible to get exact formulas and recurrences for 
the computation of quantities which are related to cycle probabilities. In particular, we 
obtain 

1. an exact formula for the appearance of a particular labeled cycle in Q^' (r) 
for (j — 2 and for a — oo (the calculation of the corresponding cycle probabil- 
ity for H(n,p) is trivial, because the edges in H(n,p) are independent of each 

1 A property P is said to be monotone if, given that it holds on a graph G, it also holds on G + e, where e 
is an edge connecting two vertices in G. 
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other). This formula immediately yields an expression for the expected number 
of Hamilton cycles in the random graph. 



2. a recurrence relation for the expected values of the elementary symmetric func- 
tions evaluated at the eigenvalues of the adjacency matrix (as a special case, 
the expected value of the determinant of the adjacency matrix) of H(n,p) and 



3. a recurrence relation for the expected values of the permanent of the adjacency 
matrix of H(n,p) and Q^h' d \r). 

These formulas can be evaluated explicitly and provide concrete information about 
random graphs with a finite number of vertices. For example, we observe that cycles 
appear earlier in GR graphs than in the ER graph. Specifically, the edge-probability 
threshold at which the expected number of Hamilton cycles crosses 1 is lower in the 
GR graph than in the ER graph. However, the difference between the two thresholds 
reduces as n increases. A similar observation can be made about the expected value 
of the permanent. The expected value of the determinant can be very different in the 
GR and ER models, indicating that for particular values of edge probabilities, the dis- 
tribution of graphs in the GR and ER models can be very different. Another interesting 
observation is that as the edge probability is varied between and 1 , the expected values 
of the determinants of the adjacency matrix can be quite large. In effect, these expected 
values provide us some useful information about the largest possible determinant of a 
symmetric 0/1 matrix. 

2 Preliminaries 

We introduce some notation and summarize some well known results to be used in the 
subsequent sections. 

We use G n to denote a random graph on n vertices (in one of the models described 
above). Then Aq u = [dij(G n )] is the adjacency matrix of G n , which is a symmetric 
random matrix with 0/1 entries (the entries of this matrix are correlated if G n is the 
GR random graph). 

Let R and C represent the sets of real and complex numbers respectively, and let 
R d , C d denote the d— dimensional spaces of real and complex d-tuples. The set of 
integers is represented by Z, and Z d is the subset of R d consisting of d— tuples of 
integers. Elements of these spaces will be denoted by bold letters such as x, y, u>. Each 
x in any of these spaces is a d— tuple (xi, X2, ■ ■ ■ xj). We will use 1 e Z d to denote the 
d— tuple with each of its entries being 1. If x = (xi, x 2 , . . . xj) and y = (j/i , J/2 • ■ • J/d) 
are two elements of these spaces, then the inner product x.y is ^ x jUj- The L a norm 
for these spaces defined in the usual way, and for x, || x \\ a denotes the L a norm of x. 
If S C R d , then S5 is the indicator function of S, so that 



QZ' d) (r). 




1 if x 6 S 
otherwise 
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For an absolutely integrable function / : R d — > R, the Fourier transform / 
R d -> C is defined as 



/M = / e— - x /(x)^(x) 

where <i/i(x) is the volume element in R d at x. Further, if /(x) = /(— x) for all 
x G R d , then /(w) = /(— to) for all u e R d , and / always takes on real values. If / 
is an absolutely integrable function with bounded support, and we define 

f p (x) = £ /(x-u) (1) 
uez d 

then f p is a well defined periodic function, that is, 

/ p (x + u) = / p (x)forallueZ d (2) 
which can be expressed by a Fourier series of the form 

/ p (x) = /( 27ru ) e2 " 4 lu ( 3 ) 

uez d 

If /, g : R d — > R are two absolutely-integrable functions, the convolution / * g is also 
absolutely-integrable and is defined as 



(/*.9)(x) 



f /(uMx-u)d/x(u) (4) 



and the fourier transform of / * g is fg. 
For r > 0, The set 

tfd,a,r(u) = {x 6 R d : || x - u || CT < r} (5) 

is termed the a— ball of radius r in R d , centered at u. The volume of Bd,<r.r(u) is 
denoted by Vd, CT ,r- Clearly, 

Vi.oo.r - (2r) d (6) 

For cr = 2 HI 

^ 2 - r = r(i + d/2) (7) 

where T is the gamma function. The surface area of B £ j IT)T .(u) is denoted by Ad, a ,r, 
and it is easy to show that Ad t oo,r = 2d(2r) d ^ 1 and that Ad^.r = dVd.2,r/r. In 
Qri' d \r), let Pd,cr,r be the probability that two vertices i, j are connected. Clearly, if 

< r < 1/2, /3 d ]J tr = V d>a , r . 

The Bessel's function of the first kind [ 15] with parameter v is denoted by J„. The 
following result is well known: 

A ( \ fn W2 J d/2{r || W || 2 ) 

V II w lb 



4 



3 The probability that a particular labeled cycle ap- 
pears in G n 

A labeled cycle in G n of length q < n is a sequence of vertices yi , y2 , • • • y q such that 
{yi,yi+i} G E(G n )fori = 1,2,... g- 1, and {y q ,yi} £ £?(G„). Let6(G„,g) 
denote the probability that this labeled cycle is present in G n . In both the GR and ER 
graph, this probability does not depend on the particular labeled cycle whose existence 
is in question. Thus, when G„ is either an ER or a GR graph, 

Q(G n ,q) =Q(G m ,q), n,m>q. (9) 

When G n — H(n,p), Q(G n ,q) can be calculated very easily. Let n > and 
1 < q < n. If G n — H(n,p), then the existence of a q— cycle in G n implies the 
presence of q edges if q > 2, and q — 1 edges if q = 2. In the ER random graph 
H(n,p), the presence of an edge is independent of the presence of the others. Thus, 

e(H(n,p),q)= | P pq *l = l (10) 

In the case of the geometric random graph Q^' d \r), things are more complicated 
because the edges are not necessarily independent. Our main result is an exact charac- 
terization of 0(Q n a ' d \r) for any a, d. 

Theorem 1 Let < r < 1/2, and q > 1. Then 

^ d \r),q) = \ ^ * (27rm) (») 

[ Lmez" "B £i ,„, r (o)^ 7^m - , 11 9 > 1 

Proof: Let xi , X2 , . . . x 9 be the g > 1 random points which form the labeled cycle of 

length q (these points are uniformly distributed in T<j). Then, Q(Qi?' d \r), q) is equal 
to the probability that for i = 1, 2, . . . q — 1, 

|| x 2 - x 4+ i < r (12) 

and || x g — xi Ho- < r. Clearly, if q — 2, then the required probability is just (3d,a,r- 
Assume that q > 2. We decompose Q(Qn' d \r), q) as follows: 

&(Q { n' d) (r),q) = Prdlxi-Xj+i \\ a < r, i = 1,2, . . .g - 1, and || xi - x, || CT < r) 
= Pr (ll x i - x 9 IU< r I II x » - x »+i 11(7 < r, i = 1, 2, . . . q - 1) 
x Pr(|| x, - x 4+1 || CT < r, i = 1, 2, . . . q - 1). (13) 

Clearly, since we are looking at i.i.d. points on the unit torus T\, the events || xi — 
x 2 ||<7< r, || x 2 — x 3 || CT < r, . . . || x<j_i — x«j || CT < r are independent of each other, 
and the probability of occurence of each is fid,<T,r- Hence, 

Pr(||a;i-a; i+ i || ff < r, i = 1,2,... = /3f^- (14) 
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Thus, we can write 

Q{Q ( °' d \r),q) = A^»x/3^ r (15) 

where 

A d M, q (r) = Pr(|| xi - x q \\ a < r / \\ Xi - x i+1 || CT < r, i = 1, 2, . . . q - 1). 

We can interpret A^o-^r) in the following manner. Consider a random walk in R d 
starting from the origin wi = 0. A point Ui is chosen uniformly in the ball Bd. a ,r(0). 
The walk then moves to W2 = wi + ui. Continuing in this manner, if the current point 
is Wk, the walk moves to Wk+i = wt + Uk where Uk is chosen uniformly in the ball 
Bd,a,r{0)- Since all points m e Z d map to the origin in the unit torus, 

Ad,a, q {r) = Pr ( w q + m e B d .a,r(0) for some m e Z d ) (16) 

Each Ui is generated uniformly from -Bd,<r,r(0), and thus, the probability density 
function of each Ui is 

P«(x) = --^ (17) 

Pd,a,r 

Then, the probability density function of Wk is the k — 1 fold convolution 

p fe (x) = (p„ *p„ * ... *p«)(x) (18) 
and the Fourier transform of pk is 

Pk(u) = j (19) 

y Pd,a,r J 

Define the periodic function s r (x) as follows 

«r(x) = E B d ^(r)(x-m) (20) 

mez d 

Then, since r < 1/2, 

Ad, a , q {r) = I s r (x)/> g (x)<fyt(x) (21) 

The periodic function s r (x) has a Fourier series representation 

s r (x) = c m e 2 ^ m - x (22) 

m£Z d 

with 

c m = S Sdi<Tr(0) (27rm) (23) 

Thus, 

A d .,, q {r) = [ V Cm p g (x)e 2 - m - x ^(x) (24) 
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Observe that p q (x) = when || x || CT > qr. Thus the integral in Eq. d24l > can be 
considered to be over a compact set, and the order of summation and integration can 
then be exchanged 1161 . and we can write 

A d ,o, q (r) = V / c mPq (x)e 2mmx dfi(x) (25) 

For any absolutely integrable / : R d — > R, we have 
f /(x)d/i(x) - /(0) 

Also, by the frequency shift property, the Fourier transform of /(x)e la x is /(w — a). 
Using these facts, we obtain 

A d ,a, q {r) = ^ S Bd ^ r( o)(27rm)p g (-27rm) (26) 

mGZ d 

From Eq. (JT9J and Eq. d26l >, the theorem follows. □ 

Using TheoremQ] we can obtain series representations for 9 in terms of the Fourier 
transform l Bg d a r (o)(w). This Fourier transform is relatively easy to compute for a = 
oo and for a = 2. 

Corollary 1 Let n > 0, < r < 1/2, and 1 < q < n. 



&(Ql?' d) (r),q) - \ { { J )dq (1 + 2 ^- i(sinc(2 ^ fcr))9) d if q >2 (27) 

Proof: Since (3d,oo,r = (2r) d , the first part of Eq. (|27| | (for q = 2) follows from 
TheoremQ] 

Assume that q > 2. Since we are using the norm, each of the d projections of 
the points xi, X2, . . . x g must induce a cycle in T±. Since the projections are indepen- 
dent of each other, it follows that 

e(Qt' d) (r),q) = (e(Q^ 1 \r),q)) d . (28) 

It is easy to see that 

a Bi lOO , r (0)M = 2rsinc(wr) (29) 

Using Eq. ( 1291 and Eq. ( 1281 together with TheoremQ]we obtain the required expression 
(we have used sinc(x) = sinc(— x) to rewrite the series). □ 

Corollary 2 Let n > 0, d > 1, ami 1 < q < n. Then 

( V d , 2 ,r q if q = 2 

e(Q "' d)(r) ' g) = { v&r + ( 2 -) d9/2 ££i Mk) { ^(0^p) q ifq>2 

(30) 

w/zere ipd{k) !S number of solutions x G Z d fo f/ie equation || x 1 1 3 = fe 
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Proof: The proof follows immediately from Eq. ([8]l and TheoremQ] □ 

Remark: In order to compute the series in Eq. (f30b , we need to evaluate the function 
ip(l(k). The following recurrence can be used: 

f 1 if k = 

^l(fc) = < 2 if k ^ and fe = m 2 for some m G Z 
[ otherwise 

and if d > 1, 

Q<m<\/k 



4 The expected number of Hamilton cycles in Qn (r) 

The Hamilton cycle problem in geometric random graphs has been studied in 1171 . in 
which the authors show that the threshold for the existence of a Hamilton cycle in a 
geometric random graph (in the unit cube) is the same as that for 2-connectivity. The 
number of Hamilton cycles in a random grapli [ 18 1 also shows a sharp thresholding 
property. 

Using 9(G„, n), we can directly get the expected number of Hamilton cycles in 
G n . Denote the expected number of Hamilton cycles in the random graph G„ by 
r(G n ). For n > 2, the number of labeled Hamilton cycles in a complete graph on n 
vertices is (n— l)!/2. It follows that, for n > 2, 

r(G n ) = Q(G n ,n) (n-l)!/2 (31) 

because the probability of each such labeled cycle being present is 0(G„, n). 

Consider the threshold for G n defined as the smallest edge-probability such that 
T~(G n ) > 1. We can use Corollaries [T] and |2] to compute this threshold when G n — 
Qn (f) an d G„ = Q^' d (r), and contrast this threshold with that for the ER graph 
H(n,p). In FigureQ] we show the thresholds obtained for H(n,p) and Qn' a \r). The 
computed threshold for the geometric random graph is lower than that for the ER graph. 
However, the difference between the two thresholds reduces as n increases. Asymptot- 
ically, the threshold for the appearance of a Hamilton cycle seems to be similar in the 
GR graph and the ER random graph (this threshold is of the order log(n)/n ifTTl ). An 
explanation for this is that as n increases, the end points of a path of length n become 
less correlated (recall the random walk argument used in the proof of TheoremQ]), and 
thus, the probability of an edge between the end points of the path is close to the edge 
probability. 

2 The random graph model used in ['18 1 starts with an empty graph on n vertices, and produces a sequence 
of graphs by adding new edges with equal probability. A threshold is then a position in the sequence at which 
a property becomes true with high probability. 
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Figure 1: Threshold for r(G„) > 1 plotted as a function of the n for the ER graph and 
for the GR graph with d = 2, a = 2 

5 The expected value of the determinant and the per- 
manent of An 

Let F Gn (x) be the matrix xl + A Gn ■ Define the two polynomials 

A Gn (a;) - det(F Gn (x))), (32) 

and 

T Gn (x) = per(F Gn (a;)). (33) 

The polynomials A Gn (x) and T Gn (x) have coefficients which are random variables. In 
particular, the coefficients in A Gn are symmetric functions of the eigenvalues of A Gn . 
Define 

A Gn (a;) = E(A Gn (x)) (34) 

and 

f Gn (ar) - E(T Gn (x)) (35) 

where the expectation of a polynomial p(x) is the polynomial p(x) whose coefficients 
are the expectations of the corresponding coefficients in p(x). 

The coefficient of x k in A Gn (x) is the expected value of the elementary symmetric 
function of degree n — k evaluated at the eigenvalues of A Gn . In particular, the constant 
term in A Gn (x) is the expected value of the determinant of A Gn , so that the expected 
value of the determinant of A Gn is A Gn (0). The coefficient of x k in f Gn (x) is the 
expected number of cycle covers across all subgraphs of G n with n — k vertices. Also, 
the expected value of the permanent of G n is f Gn (0). 
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ER Random Graph + 
GR with d=2, delta = 2 x 



There is a strong connection between cycles and permutations, and between per- 
mutations and determinants (and permanents). We expect that the characterization of 
6(G„, q) will help determine the behaviour of the determinant (and permanent). More 
concretely, we show that 

Theorem 2 Let G n be a random graph on n > vertices ( G n is either the ER graph 
or the GR graph). Then, for n > 1, the polynomials A Gn {x) and f Gn (x) satisfy the 
recurrence relations 

A Gn (x) = xAg^Jx) + f^i-iy- 1 — - Q(G n ,q)A Gn _ q (x) (36) 

and 

n _ II 

f Gn (x) = xTg^x) + V ^— Q{G n ,q) t Gn _ q {x) (37) 

n - q - 

with initial conditions A Go (x) = f Go (x) = 1. 

Proof: We start with the following formulas for the determinant and the permanent. If 

B = [bij] is an n x n matrix, then 

n 

det(B) = ^(-ir 9 " (CT) n & -« 08) 

<tS5„ i=l 

and 

n 

per(S) = J2 U b ^) ^ 39 > 

o-GS„ i=l 

where S n is the group of permutations of {1, 2, ... n}. 

Each permutation a e S n can be uniquely decomposed into a set of disjoint cycles 
on {1, 2, . . . n}. Each cycle C in the disjoint cycle-decomposition of a permutation is 
of the form (iii2 ■ ■ ■ i q ), where a(i r ) = i r +i, r = 1, 2, . . . q — 1 and <r(i q ) = i\. The 
sign of the cycle C is sign(C) = (— l)' c ' _1 , where \C\ is the number of elements in 
C. The sign of the permutation is then the product of signs of the cycles into which a 
is decomposed. We will say that the pair £ C if i,j are consecutive elements in 
the cycle C (io is considered to be after i q ). Then, given a, we have 

n 

i[b Mt) = n n b v < 4 °) 

For a cycle C, We define 

(ij)ec 

Then, 

det(B) = E II (-1) |C|_1 ^b(C) (42) 

cr£S„ C6<t 
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and 

per(B) = ]T J] w b{C). (43) 

treS n CGcr 

Let B — Fc n (x). For a cycle C = . . .i q ) in some permutation, we see that if 
q > 1, then 

E{w B {C)) = Q(G n ,q) (44) 

and if q = 1, then 

E(w B (C)) = x. (45) 

For convenience, we set 0(G„, 1) = x. 

Also, if C\ 1 C*2, . . . Ct are vertex-disjoint cycles in G n , then the presence of Cj is 
independent of the presence of Cj for j ^ i, and 

t t 
E(H w B {Ci)) = H E(w B (Ci)). (46) 



It follows that 



Similarly, 



a g„(x) = e n(- i ) |cHie ( G -i c i)- < 47 > 



r Gn (x) = E II e ( G «-l c l) («) 

The counting of permutations a 6 S* rl can be carried out by fixing the cycle C which 
contains 1 and counting permutations of elements not in C. For 1 < q < n, Let D q be 
the set of cycles of length q which contain 1. We observe that 



\D q \ = 



n- 1 

q-i 



because each cycle in D q is determined by the choice of q — 1 elements (other than 1) 
out of n — 1 elements, and there are q — l! distinct cycles on q elements. 

Let N = {1, 2, . . . , n} and let P(A) be the set of permutations of the set AcN. 
Then, we can write 



E E e(Gn, E II (" 1 ) |D| " 1 e(G B , (50) 

g=l \C£D, \<reP(N-C) B€tr / / 

where the innermost summation over P(A) is taken to be 1 if A — </>. Since C = g 
for each C € Z? 9 , we can rewrite Eq. (l50l > (using Eq. (O to replace 0(G rl , Z?|) by 
9(G n _ g , |£>|)) as 

E( n ~\ ) (q-iy.i-iy-'eiG^q) I e nt- 1 )' 1 ^ 10 ^"-^!) 

<?=i^ ' \o-eP(N-c) -Deo- 

(51) 
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The inner summation in Eq. ( BTT i is just A(? n (x), and thus, the recurrence relation for 
A.G n (x) follows. The recurrence relation for Ta n (x) can be shown to hold in a similar 
manner, completing the proof of Theorem[2] □ 

Remark: The result in Theorem|2]holds for any random graph G n in which the proba- 
bility of appearance of a labeled cycle depends only on its length and the probability of 
appearance of a set of vertex-disjoint cycles is the product of probabilities of appear- 
ance of the elements in this set. 

For n > 0, < k < n, let F n k(ti, . . . t n ) denote the elementary symmetric 
function 

Fn,k(t±, t2, ■ ■ ■ t n ) = ti l ti 2 ...ti k (52) 

{il,i2,...ik}e{l, 2, ...,«} 

For k = 0, define F n ^ = 1, and define F n ^ — if n < k or if k < 0. Now, let F n ^ 
denote the expected value of F n k evaluated on the n eigenvalues of Aa n . Then, the 
expected value of the determinant of A& n is just F n>n , Then, we have the following 
corollary of Theorem|2] 

Corollary 3 For the random graph G n , ifn>0, and < k < n, then 

n n — l! 

F„,fe = F n _ hk + V (-l)' 1 e(G n ,q) F n _ q>k _ q (53) 

Proof: Follows from Theorem|2]by noting that the coefficient of x k in A is F n> n-k- O 



Note that in both models, if the edge probability is 1, then 0(G„, q) = 1, and G n 
is always the complete graph, so that the expected value of the determinant of G n is 
(— x (n — 1). Using Theorem[2] we obtain the following identity for n > 0: 

n = 1 + ~ x (( n " - 1) (54) 



n — o! 

g=2 y 



Also, the permanent of the complete graph on n vertices is the number of derangements 
of the set N = {1, 2, ... n}. Thus, the recurrence proved in Theorem [2] yields the 
following identity for the number of derangements d n of N 



E 



77-1! 

dyi—qi (55) 



q =2 71 - 



with the initial conditions d\ = 0, and do = 1. 

We use these ecurrence relations to compute these expected values for n < 20 
in the GR and ER modelfl Some interesting conclusions can be drawn from these 
calculations. 



3 The recurrence relations were directly computed using long double precision arithmetic. For higher 
values of n one would need to use higher precision arithmetic. 
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Figure 2: The expected value of the determinant plotted as a function of the edge 
probability for n = 20 in the ER and GR (with d = 3, a = oo) models. 

Consider the plot in Figure [2] in which we compare the behaviour of the determi- 
nant of G20 as a function of the edge probability. The graph has been plotted for the 
ER graph and for the GR graph with d = 3. The behaviour of the determinant in the 
two models is quite different, and clearly, so is the distribution of G n . 

In Figure [5] we show a plot of the expected value of the permanent of Aq„ (for 
n = 20) as a function of the edge probability in the ER and GR (d = 1) models. We 
can also define a threshold for the expected permanent as the smallest edge probability 
for which the expected value of the permanent is > 1. A comparison of this threshold 
for the GR and ER graphs shows that this threshold is lower for the GR graph, but 
the two thresholds come closer as n increases (see Figure 0J. Thus, the permanent of 
the GR graph grows more rapidly than that of the ER graph. This is expected since a 
labeled cycle is more likely in the GR graph. 

5.1 Graphs with large determinants 

Looking at Figure [2] we see that for intermediate values of the edge probability, large 
magnitudes appear in the plots of the expected value of the determinant. For instance, 
we observe that, in the ER random graph with n = 20, the largest absolute value of the 
determinant is 3787.81, and this provides a lower bound on the maximal determinant 
of a symmetric 20 x 20 0/1 matrix . 

For a general (non-symmetric) n X n 0/1 matrix, the determinant is bounded above 
by (?i+l)(" +1 )/ 2 /2" [19]. The number of (possibly non-symmetric) rixn 0/1 matrices 
which achieve this bound is also known for n < 9 [20|. However, similar characteriza- 
tions of the determinants of symmetric 0/1 matrices are not so common. For example, 
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Figure 3: The expected value of the permanent plotted as a function of the edge prob- 
ability for n = 20 in the ER and GR (with d = 1, a = oo) models. 



ER Random Graph + 
GR with d=2 : delta = infinity x 



12 14 16 18 



Figure 4: Threshold for the expected value of the permanent plotted as a function of 
the n for the ER graph and for the GR graph with d = 2, a = oo 
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in ||2T1 . the authors show that for n > 7, the maximal determinant of the adjacency 
matrix of a (n — 3)-regular graph on n vertices is {11 — 3)3^"/ 4 l -1 . For n = 20, this 
works out to be 1377 which is less than the largest observed determinant value in the 
evolution of H(20,p). 

Thus, the recurrence formula for the expected value of the determinant seems to 
provide some useful information about the maximal determinant of a class of sym- 
metric 0/1 matrices (in effect, we have a lower bound on the largest value of such 
determinants). Also, if the expected determinant is large, then it may be possible to 
find a symmetric 0/1 matrix with large determinant by using a Monte Carlo sampling 
approach. An estimate of the second moment of the determinant of the random graph 
will throw more light on this possibility. 

6 Conclusions 

We have derived an exact characterization of the probability of existence of a labeled 
cycle in geometric random graphs on a unit torus with an arbitrary number of dimen- 
sions, and with an arbitrary L a distance metric). This cycle probability can be calcu- 
lated in terms of the Fourier transform of the indicator function of a ball in L a . Explicit 
expressions for this Fourier transform can be easily computed in the a = 00 and a = 2 
case. 

From the cycle probability, one gets the expected number of Hamilton cycles in the 
geometric random graph. These exact expressions complement the asymptotic thresh- 
old results for the existence of Hamilton cycles in geometric random graphs (as in ifTTl ). 
We observe that as the edge probability increases, a Hamilton cycle appears earlier in 
the GR graph than in the ER graph. 

The cycle probabilities can also be used to find the expected values of the determi- 
nant (and more generally, the expected values of the elementary symmetric functions 
evaluated at the eigenvalues of the adjacency matrix) and the permanent of the adja- 
cency matrix of the random graph. We obtain recurrence relations for these quantities 
and illustrate them by a few calculations. In particular, the determinant exhibits very 
different behaviour in the two models. Also, large magnitudes of the determinant are 
observed in the evolution of the random graphs. This throws some light on the as yet 
unresolved question of the maximal determinant of symmetric 0/1 matrices. 
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